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Abstract: Consider the classical XY model in a weak random external field pointing along the Y 
axis with strength e. We study the behavior of this model as the range of the interaction is varied. 
We prove that in any dimension d > 2 and for all e sufficiently small, there is a range L = L(e) so 
that whenever the inverse temperature /3 is larger than some /3(e), there is strong residual ordering 
along the X direction. 



1. Introduction 

In this paper we study an interesting phenomenon which has received little attention in the 
mathematical physics literature: random field induced ordering. The paradigmatic example is 
given by the following classical XY spin system. Let {a z } zeIi d = {a z (uj)} zeIi d = be a family 
of i.i.d. {±1} valued random variables taking each value with equal probability. Configurations 
are denoted by uj G £1 = {0, 1} Z with (fi, B, P) the corresponding probability space. Let e > 
be fixed. The (random) Hamiltonian of the model of interest is 

JT(ct) = JT<» = - a * ■ a *' ~ e £ ot z (Lj)z 2 ■ a z (1.1) 

(z,z') zeAjv 

where the sum is over the set of nearest neighbor bonds in Z, d . Here for each z G a z G S 1 and 
we take as a priori measure du(a z ) Haar measure on S 1 . Finally 61 , 62 denote the two vectors 
(1, 0), (0, 1) respectively. 

Superficially, (1.1) is at the intersection of a number of notable d = 2 phenomena: On the one 
hand, if e = 0, we have a pure XY model, which does not order (magnetically) at any temperature 
as a consequence of the Mermin-Wagner theorem [6]. On the other, if the spin space is take to 
be {—£2, £2} rather than S 1 , then the model is the random field Ising model (RFIM), which 
does not order either as explained heuristically by Imry-Ma [4] and demonstrated rigorously by 
Aizenman-Wehr [3]. 

The question then is what happens when the two models are "combined"? In the (physics) liter- 
ature it is expected, surprisingly at first, that there is in fact magnetic ordering at low temperature 
for any d > 2 and, at present, no truly rigorous mathematical results exist (in any dimension). 
The d = 2 case is particularly subtle for a number of reasons. An even more surprising feature of 
this model is the direction of the ordering, expected to occur along the Si direction. 
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Interest in systems of this type have seen a flurry of activity in the physics literature in the last 
few years; a very recent review is [8] which focuses on these effects in the quantum context, see 
also [2,9, 10] for work on classical systems. A related mechanism, not covered in [8], appears 
in [1]. There, the interesting point is that the magnetic field is determined by a random gauge 
potential. As a result the average net field strength in a cube of side length L scales as \fh in two 
dimensions as opposed to the scaling of L in (1.1). 

Since we have so far been unable to treat directly the nearest neighbor model, to gain under- 
standing of the mechanism behind the ordering we turned to its mean field theory. Some light 
is shed in this simplification as we find that the free energy landscape of the model has a double 
well, with minima at (roughly) ±p((3)ei for some positive p{f3) tending to \/l — e 2 as /3 — > oo. 
This observation is the starting point for analysis of the system (2.1) when we take lattice models 
with Kac interaction of range L = Poly(^). The result is formulated as Theorem 2.1. No attempt 
has been made to optimize the degree of this polynomial, although the discussion which follows 

_ 2 

indicates that a lower bound of L > e a is necessary to approximate lattice systems by mean 
field theories. 

We should note that the mean field theory was previously addressed in [2, 9] in somewhat 
different ways. Those treatments are not sufficient for our purposes and moreover, the authors 
focus on critical behavior in the (e, j3) plane. As we will argue, the link between this and the 
behavior of lattice systems seems tenuous, at least in the weak field regime. 

To provide context to the free energy functional will we introduce in a moment consider the 
formulation of (1.1) on the complete graph on N vertices: the Hamiltonian becomes 

3?{<y) = -— ^2 &z ■ (?z' ~ e ^2 a> z e 2 ■ &z- (1-2) 
z,z'e[N] ze[N] 

It is not difficult to compute <j>, the negative of the large deviation rate function for the vector 
observable (M^,M^), where denotes the spatial average spins a z over {z : a z = ±1} 
respectively. For m € M 2 , let S(m) = inf heR 2 (G(h) — m ■ h) with G(h) = log J §1 Av(a)e a ' h 
denoting the log moment generating function of v. Let 

<fi(m + ,m~) = 4>/3 ie (m + ,m~) := --||m||| - |e 2 ■ (m + - m~) - ^( S ( m+ ) + s ( m ~)) 

where fh = \{m + + m~) and f(m + ,m~) = 4>(m + , m~) — inf( m + m -) 4>(m + , m~). We use 
the notation 

\\(m + ,m~) - (m^,rriQ )|| = max(||m + - m^|| 2 , \\m~ - tUq || 2 ) 
with || • || 2 denoting the usual Euclidean norm. Then 

where is the Gibbs state associated to (1.2) and convergence occurs for all (m + , m~) w-a.s. 
The main properties of 4> are summarized in the following theorem. 

Theorem 1.1 (Low Temperature Mean Field Phase Diagram) There exists eo > and (3q > so 

that for all e < €q, (3 > (3q, the free energy functional f(m + , m~) has precisely two minimizers 
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±(m . m ). These minimizers are characterized by 

||m + || 2 = ||m _ || 2 = p, 

m + • e\ = mT ■ e\ = cos(0). 
where, p < 1 and 8 G [ir/2,ir/2] satisfy the mean field equations 

sin(0) = -, 
P 

p = —d p S{pe\) ( the maximal solution, which is nonzero). 

In particular, \p\ is bounded away from Ofor all (3 sufficiently large and e small and consequently 
8 = 0(e). 

Further, we have the following stability estimate: 

4>(m + , m~) — 0(m + , m~) > 

c(e ,/3o)||(m + ,m~) + (m + ,m~)|| A \\(m + ,m') - (m + ,m~)|| A e 2 (1.3) 

with all (non-minimizing) stationary points (m,Q , ) satisfying 

e 2 

0(m+,m o ) - </>(m + ,m ) > — . 
The above theorem says that free energy <p has exactly two minimizers and that the height of 

2 

the minimax barrier between the two minimizers is of order y uniformly in j3 large. It also says 
that these two minimizers are transverse to the direction of the randomness so that fluctuations of 
the local fields do not favor one of these minima over the other. This gives a stability not present 
in the RFIM. 

To make the mechanism for ordering even more transparent, we may attempt to interpolate 
between the random field transverse field XY model and the RFIM by considering Z n Clock 
Models. The spin space is the set of nth roots of unity on the unit circle, viewed as vectors in 
M 2 , with Hamiltonian as in (1.1). If n = 4, the randomness forces the system to behave much 
as in the case of the RFIM, as can be seen purely from the consideration of ground states. One 
may then wonder what happens as n is increased. We find, at least for the ground states, that for 
each e > 0, there is a crossover behavior from RFIM to the random transverse field XY model 
occurring at n ~ K That is, if n « ^ then ground states oscillate between ±pti2 where as if 
n » 2, there are a finite number ground states, stable under the noise, all roughly parallel to Si. 
In particular, in contrast to statements made in [9], a spin space with truly continuous symmetry is 
not necessary for the effect of interest to occur. All that matters is that the spin space has enough 
freedom to take advantage of local fluctuations. 

We would also like to point out that one can perform a mean field analysis in case P(a x = 
1) = p / \ and a similar picture emerges in the mean field theory. In particular letting q = 1—p, 
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the mean field equations become 



pm ■ ei 
pm + ■ &2 — • &2 



qm ■ ei 



€ 



1 



d p S(p+ei) 



(maximal solution) 



P+ 



2pp 



1 



d p S(p-ei) 



(maximal solution) 



P- 



2qP 



as can be seen by following the proof of Theorem (1.1). These equations imply that at low 
temperatures, the picture presented above in the unbiased case persists-there are two minimizers, 
symmetric with respect to the F-axis. The location of the average of the two components of a 
minimizer has a non zero Y component with sign and magnitude determined by the bias. 

Let us return to requirement L = Poly(-). Fluctuations of the local fields are of order y/N 
typically, and can change the finite volume low temperature free energy landscape by O(-^U). 
Thus the most likely order parameter values at finite volume can have 0-free energy which differs 
from the absolute minimum of <p by an amount of order Large deviations implies that the 

height barrier in a system of N vertices is of order e 2 N, so the effects of field fluctuations on 
macroscopic observables will only be suppressed if Ne 2 >> eVN. This leads us to a fundamen- 
tal requirement for the validity of the 0-mean field picture even to the model on the complete graph: 
y/N » \. Translated to Kac interactions, this means we must take the range of our interaction 
L so that L d/2 » \. 

The above discussion provides a first noteworthy point of our work. Though it is not believed 
that randomness induced ordering depends on taking a long range interaction, our work empha- 
sizes the difficulties in drawing conclusion about short range models (even in high dimension) by 
extrapolating from mean field analysis, particularly when the random field strength is weak. This 
point seems to have been overlooked in the literature. 

A second main point (somewhat counter to the first) is that we are able to give a rigorous 
example of a lattice system where the combination of a random field acting in one direction 
with a coupling that has continuous symmetry disrupts the behavior of the two d = 2 systems 
discussed above (Mermin-Wagner and RFIM). Indeed, it is known that neither of those models 
order at low temperature even when the interaction range is spread out as Poly(^). 

The rest of the paper is organized as follows. In the next section, we precisely formulate our 
main result Theorem 2.1: the existence of residual magnetization along the Si direction for Kac 
models with sufficiently long range interaction. Section 3 states the main Lemmata needed for 
the proof of the Theorem 2. 1 and on their basis provides a proof of the Theorem. Section 4 is 
devoted to the proof of Theorem 1 . 1 and other technical estimates needed regarding the mean field 
theory. There after, the paper is devoted to a justification of the Lemmata appearing in Section 3. 
The techniques used are mostly taken from the book [7], with modifications needed to treat the 
randomness and the fact that we are working with a continuous spin space. This means making 
appropriate definition of course grained contours and comparing contour energies to a certain 
free energy functional evaluated on deterministic magnetization profiles which are defined on the 
support of each contour. 
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2. Main Result 

Let S = (S 1 ) zd be endowed with the product topology and associated Borel cr-field Bq. For 
any a G S and any set A C 7h d let a a denote the restriction of a to A and Aua = du(aA) 
denote the corresponding restriction of v. Sa, <6o,a will denote the associated space of spin 
configurations and sigma fields respectively. 

Let Ajy = {z G Z d : ||z||oo < N} where ||z||oo denotes the length in R d . and introduce a 
length scale L which represents the range of the interaction. So Jl{z, z') := clJ(\\z — z'^/L) 
where J is C c °°({x G R d : \\x\\ 2 < 1}), c L x = / dxJ(||x||/L) is a normalizing constant. Let 
e > be fixed. Given cjac g S\c , the random Hamiltonian for our model is 

^n(<ta n Wa%) = ~ ^2 Jl{z,z!)o z - a z > - a z e 2 • a z - ^ J L (z, z')a z ■ oy 

z,z'eAjv zeAjv zeA N ,z'£A c N 

Notice the crucial property that, even in the presence of the random fields, the energy ^ai{pa n \va% 
is invariant with respect to simultaneous reflection of all spins about the F-axis. 

The Gibbs-Boltzman probability distribution associated with this Hamiltonian is defined by 
Radon-Nikodym derivative relative to va n '■ 

d/% n (cta n ) =d/i N N ((JA N )cKe N > a n>& v (oa n ) (2.1) 

with constant of proportionality Z^(j3, &a c n ) the (random) partition function of the system. 
Horizontal Boundary Conditions: For convenience, let be the random Gibbs state with 
— > denoting boundary conditions given by setting all boundary spins equal m := ^ (m + + m ) . 
Even though this choice is not strictly in the spin space, the model still makes sense. Further, this 
external configuration can be effectively produced mesoscopically by taking spins to be of the 
form aei ± 6e 2 where the sign varies according to the parity of the underlying lattice site. In this 
way one obtains block average magnetizations aei + 0(k~ 2 ) where k is the size of the box. Since 
we are considering Kac interactions, the effect is the same as our Horizontal Boundary conditions 
to within an error which plays no significant role. 

Given A > 0, let l> = L 1+x . Let B e z > = {x G Z d : \\z - x||oo < £>}■ A subset A of Z d will 
be said to be £>- measurable if A is a union of blocks Br> so that r G {2£> + l)Z d . Let 

Bf*> = {x G B e > : a x = ±1} 



We define the block average magnetizations 

\B 



M z = 



1 



xeB 



I p>£> I 

1 1 xeB l > 



Theorem 2.1 (Main Theorem) Let d > 2 be fixed. We can find £o, eo > and A > so that for 
every < e < eo and < £ < £o> there exists f3o(e, £), Lo(e, £, A) > Ofor which the following 
holds: If '(3 > /3q,L > Lq, then for almost every oj G ft, there exists an £ > = £ > (L)-measurable 
subset V u C Z d and an Nq{oj) G N such that: 
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(1) 

\VH A N \ < q\A N \ 

for all N> N (u). 
(2) For each z e A N with B e > (z)DV = 0, 

and 

\\(M z )^-rh\\ 2 <^ 

Explicitly, we require 

e 2 A A L~ Xd > C{L^ x -^ +2Xd + L 2Xd e-<* 2LX ), 

L~ x < g and £ > cL^^ l ^ x \ where c, C are universal constants. 

Remark 2.2 Let us note that the proof we employ below will work, with the appropriate modifi- 
cations, in the case of bias P(a 2 = 1) = p ^ \. 

Notational Convention Below the constants c, C, C\ > will always be universal in the sense 
that they only depend on d and the sup norm of V J but not on £, e, (3, L etc. Their values may 
(will) change from line to line. 



3. Course-Graining and Contours 

We define £ < = [L 1 ~ A J and = |~L 1+A ] to be two scales respectively slightly smaller and 
slightly larger than L. We shall assume 2£ < + 1 divides L and L divides 2£> + 1. 

The scales and introduce a filtration of blocks in M. d and, by taking intersections, in Z d . 
We shall say that a block B r = {x G R d : \\x — r||oo < L'} is measurable relative to the scale V 
if r G (2L' + l)Z d . For A C Z d finite, let 

iVf = {a x : x G A, a x = ±1}. 

It is standard that 

P(|7V±| - |A|/2 > A^/\A\/2) < 2e{-A 2 /A). 

where | A| denote the cardinality of a finite subset of Z d . For our purposes, it will suffice to take 
A = |A| K for some k G (0, 1/2). If A = B r C 1, d for some £<-measurable block B r , we will 



, L1UUU1V1J V ' 

We define 



use the shorter notation N r 



where 2? r denotes the £<-measurable block containing z. Note that this depends on the realization 
of the randomness. Also let 

Z < X&B r 

which does not. For D C^ d Borel measurable, we shall use the notation 

L°°(D) = {m: D ^ R 2 s.t. m is Borel measurable and esssup^ g £> llm^l^ < oo} 
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where \\m z H2 is the usual Euclidean norm in R 2 (more generally this notation is used for 2-norms 
in all Euclidean spaces R d ). We will also use ||m z ||oo = max^^ \m z (i)\ for vectors m z G M. d . 
The induced norm on L°°(D) is denoted by 



MIl~(d) = esssup^ 6D \\m z \\2- 



For any pair (m + ,m ) G L 
z G B r with B r £< -measurable 



00 (TU>d\ 



L°°(]R d ), we denote m = ^(m + + m ) and for any 



rn 



e<,± 



< JB 



dy m± 



and 



m< = -(m z <,+ + m<< 



Below, we refer to magnetization profiles (m + , m~) which are piecewise constant over the set of 
all B r which are i < -measurable (and in the domain of definition of (m + , m~)) as I ^piecewise 
constant. 

Fix A C Z d which is -measurable. Given a spin configuration a G S, and for any z G 
A,zeB r such that £? r is £<-measurable, we introduce the phase variables: 



7/ = 7/f 



fl if||(^<' + ,af<-)-(m+,m-)|| 00 <e, 
-1 if||(^<' + ,^<'-) + (m+,m-)|| 00 <C, 
otherwise. 

• For any z G A, z G B r with £? r -measurable 

1 if 7/y = 1 for all y G B r , 
—1 if if rjy = — 1 for all y G -B r , 
otherwise. 



6 = 9*1 = 



These phase variables are extended to the set of boxes B r C A c with dist(£? r , A) < ^> 
via 



v = vl = vl z 



1 ifllo - . 



mi 



-1 if Her/ + m||oo < £, 
otherwise. 



and similarly for 6. We will always work under conditions in which 9 = ±1 over con- 
nected components of A c . 
For any z G A, z G -B r with £? r ^> -measurable 



e = ef = e« = 



1 if 9y = 1 for all 7/ : ||z — j/H^ < 
— 1 if 9y = 1 for all y : ||z — j/H^ < 
otherwise. 



8 



N. CRAWFORD 



We emphasize that depends on A and the boundary conditions. Another way of defining 6 is 
to say that an ^> -measurable block has nonzero value of 6 if that block and all neighbors (in the 
£°° sense), including boundary boxes, have same nonzero value. 

For any -measurable region D in M. d , we may extend all of these notions to deterministic 
magnetization profiles (m+,mj) G L°°(D) x L°°(D). Note that the definition of requires 
the auxiliary input of a fixed boundary condition fho iZ G L°°(D C ). The boundary condition used 
will be clear from context. 

For any set A C Z d we can associate a subset A C M. d which is the union of boxes of side 
length 1 centered at the elements of A. We shall say that A is connected if A is (note this is NOT 
the same as connectivity in A c decomposes into one infinite connected component Ext(A) 
and a number of finite connected components (Intj(^4))^ 1 with Int(A) = U™iIntj(A). Let us 
denote the enlargement of a set A C M d by 

5(A) = U{B r e > -memmHe:dist(Br,A)<e<}^ r 

where dist(I? r , A) is the Hausdorff Distance between sets in R d in the metric. From here 
we may introduce 5(A), lnti(A), etc. by taking intersection of each defined set with 'L d . The 
closure of a set A C 7L d is defined to be c(A) = 5(A) U Int(A). It is standard that the connected 
components of R + = {z G 7L d : 6f = 1} and Rr = {z G Z d : @l = -1} are separated by 
connected subsets of R° = {z G Z d : @l = 0}. 

Definition 3.1 A contour F is defined to be the pair (sp(r), 6r) where sp(T) C Z d is connected, 
iy-measurable and 0r(z) is an £ < -measurable { — 1, 0, l}-valued function on sp(T) which gives 
the values of the phase specification on T. In the previous definitions, whenever the set A in 
question happens to be sp(T) we will write 5(F), c(T), etc. Let N r = |<5(r)|/(2f> + l) d , so that 
N r G N. 

Let us introduce the notation 5 ext (T) = 5(T) n Ext(A) and <5f n (T) = 5(T) n Inti(A). These 
are evidently disjoint. By definition of T, each of these sets is connected (in our sense) and 
disconnected from the rest. 

We shall denote by 

X(r) = {a G S : T is a contour for a} 
:= {a : sp(T) is a maximal connected subset of R°(a) and z (a) = 9r(z) on sp(T)}. (3.1) 

We shall say that F is a contour for a if a G X(T). 

By our definitions, specifying that F is a contour of a lets us recover the values of Q z (a), 6 Z (a) 
on 5(F) (see [7] for details). This convenient property allows us analyze systems of contours 
without worrying about the microscopic spin configuration far away from the contour. Two con- 
tours Ti, F2 are said to be compatible if <5(Ti) n sp(T2) = and 6r 1 = 9r 2 on tne domain of 
intersection of 5(Ti), 6^2). 

So far we have considered contours at the level of spin configurations. We would like to be 
able to show that under certain finite volume Gibbs measures, a contour costs e ~ c ^' e< ) Nr , the 
constant c being made large by appropriate choice of /3, £, L. In general, such an estimate will 
NEVER be true uniformly in the presence of randomness, but we can hope that for the family of 
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-measurable blocks {B r }, large fluctuations of the variables — \ B r \/2 are sufficiently rare 
so that we can still extract e~ c '^' e< ^ Nr in cost from MOST contours. 

Thus we turn to the interplay between spin configurations and randomness. We introduce 
(more) phase variables associated to the randomness - wed - restricted to boxes B r C 7L d with 
side length 2£< + 1. Let 1/2 > k > be fixed. For B r ^-measurable and z G B r define 



and for any z G B r , with B r l>- measurable, 



1 if \N± r - \B r \/2\ < |S r |V 2 +« 
otherwise. 



1 if <£> y = 1 for all y £ B r 
otherwise. 



Definition 3.2 Let us say that an t-> -measurable set S is ( n,p)-clean, or just clean, if 

S clean :={zeS: ~ e> , z = 1} 

satisfies |S' clean |/|S'| > 1 — p. Otherwise, S is called dirty. A contour F will be called clean if 
6(F) is clean. 

Let 

A = {Y C Z d : Y is ^-measurable, connected, and (n,p) dirty} 

Let V := Uy e x:c(y). Note that these definitions only depend on the realization oj and not on 
possible spin configurations. 

Given a contour F, F* will denote the (f2, £>, P) event 

T* = {6(F) is clean and c(F) is not strictly contained in V} (3.2) 

If T* occurs, T will be called a *-clean contour. Given a spin configuration (a\ N , <t\ c n ), let 

x(r*, . . . , r^, r m+ i, . . . , r m+n ) = riiX(rj) 

where (Ti, . . . , F m ) satisfy the event defined in (3.2) and (T m+ i, . . . , F m+n ) do not. Otherwise 
we define the right hand side to be the empty set. As a variation of standard definitions, let us say 
that 

(r*, . . . , rj^, r m+ i, . . . , F m+n , u) 

are *-compatible if 

x(r*, . . . , r m , r m+ i, . . . , F m+n ) / 0. 

Lemma 3.3 There exist 5, eo > so that for any p, X,k G (0, |), < £ < 5 so that L~ x < | 
and £ > cL - ^ 1- ^) log L and if e < eo, there exists f3 e so that if f} > j3 e then the following holds: 
Let N be fixed and consider the event 1L(F\ , . . . , F* m , F rn+ \ , . . . , F m+n ) with sp(Tj) C An- Then 

^(X(r;, . . . , v* m , r m+1 , . . . , F m+n )) < e-'^i ^ 
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where 

q = C 1 pL d ^{e 2 Ae^L- Xd - 

C ( L \-l+2\d + epL 2\d + eL (-d/2+dK)(l-X)+2\d + L 2\d e -ce 2 L* + ^-1^ 2Xd- § logL ) J 

(3-3) 

wifft c, Ci, C > universal constants. 

This lemma is proved in several parts below. Theorem 2.1 is then completed via the following 
Peierls contour counting argument. 

Proof of Theorem 2.1. Fix x G Z d . Throughout this proof, let B(x) denote the f>-measurable 
block containing x. Consider the collection of bounded ^-measurable connected subsets of Z d 
containing x: {Y}y 3x . For each such Y, we first estimate the event that Y is -dirty. We 
have: 

F(Y is (/e,p)-dirty) < 2 NY e~< H dpNv . 

where Ny is the number of £> -measurable blocks in Y. Now the number of -measurable 
connected sets Y containing x with Ny = r is well known to have the asymptotic Oq for some 
fixed, dimension dependent constant cto. Thus 

F(B(x) is in some (K,p)-dirty Y) < ^(2a ) r e _c '< <ipr . 

r>l 

Modifying the estimate slightly, via the discrete isoperimetric inequality 

F{B(x) is in c(Y) for some (k, p)-dirty y ) < C ^ r^^ 1 ) (2a ) r e- rf < KV . (3.4) 

r>l 

where C is a universal constant coming from the isoperimetric bound. 

Next we need a correlation bound. Let A(x) = {B(x) is in c(Y) for some (K,p)-dirty Y}. 
Using the fact that the events {c(Yi) is (n,p) — dirty} are independent if c(Y\) n c(Y 2 ) = we 
have 

P(A(a;i), A(a; 2 )) - F(A( Xl ))F(A(x 2 )) < 

P(B(xi), B(x 2 ) are in c(Yi), c(Y 2 ) for some (/c,p)-dirty YI, Y 2 with c(Yi) n c(Y 2 ) / 0). 

(3-5) 

Then, estimating as in (3.4), the right hand side is bounded by C£ d > e~ ci '< dpdiste >( x ' y " > where 
dist^ > (x, y) denotes the minimal number of blocks in an -measurable block path from B{x) 
toB(y). 

These bounds imply a constraint on parameters: namely 

L 2Kd ^-^p>c\oga 

We will assume from now on that n= \,p = L~»^ 1 ~ A - > and A < |, so that the inequality holds 
for all L > Lq for some Lq G N. 
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Recall that 

A = {Y C Z d : Y is (n,p) dirty, ^-measurable and connected}, 

V = U Y eAc(Y). 

Letting V N = Vf] A N , V N = U xe \ N A(x) so E[|Djv|] < C\A N \e~ c£2 < d P and, via (3.5), 

Var(\V N \) < C^/\A~M 

By taking N = 2 k applying Chebyshev's inequality to estimate deviations of |Pjv| and then the 
Borel-Cantelli lemma along this subsequence, we have that, for almost every oj € 0, there is 
N (u) e N so that 

^4 < C£ d >e - ce < dp < CL 2d e- cLl 
\^n\ 

for all N > Nq(uj). We have used here that T>n increases with N. 

The rest of the proof is an application of Lemma 3.3. As p = ( X ~ A ) we may now assume 
the dominant contributions to qL~^-~^ d are 

/3(-e 2 A e A L~ Xd + L ¥ x -^ +2Xd + L 2Xd e- ce2LX ) + L 2Xd . 

We thus require 

e 2 A e A L~ A<i > 2(L§( A - 1 )+ 2Ad + L 2Xd e~ ce2 ^ X ) 

which implies g > ciL d ( 1_A )/3e 2 A£ 2 f\L~ Xd - c 2 L 2Xd . Clearly for each d we may find A = X(d) 
and then L = L(X, e, £) so that this holds for all L > L(X, e, £). Then we see that g can be made 
arbitrarily large by the appropriate choice of /3 with the rest of the parameters fixed. 
Fix x € Aat so that .Er> (x) C\V = and consider the event 

which is a subset of Sa n - By definition of 6, and the — ► boundary condition, there exists a 
largest contour T so that B > (x) C c(T). Moreover, since I? £> (i)nD = 0,r must be (n,p)- 
clean. Decomposing {Ob 1 >{x) / 1} into disjoint subsets according to this largest contour we 
have: 

r,w *-compatible: _B*> (x)Cc(r) 

By Lemma 3.3, 

r>l 

as long as g > 2 log(2ao). The theorem now follows easily. □ 
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Proof of Theorem 1.1. Let us use the polar coordinate parametrization = (p±, 9^). Although 
there are cleaner ways of deriving what we need, as kindly pointed out by D. Ioffe, we will stick 
with direct computations in these coordinates. 

Evidently all minimizers have the property that either 9 ± 6 [— ir/2, ir/2] or 9 ± G [vr/2, 3ir/2]. 
By symmetry, it is enough to treat the case 9 ± G [— tt/2, tt/2]. We may also assume 6 + > 0~ . 

In polar coordinates 

= ttPi,^) = ~\{p\ + P-+ 1P+P- cos(# + - n) - € 2 (p + sm(9 + ) - p_ sm(9-)) 

-^(5(p+ei) + 5(Mi) 
= l(9f}(p+) + 9(i(p-)) ~ \p+P- cos(6 + - 9-)) - € -{p + sin(0+) - p_ sin(O) (4.1) 

where 

Considering the angular gradient, at any stationary point we have that 

= V e+ <P = sin(#+ -9-)- 6 -p+ cos(#+) 

= V e -<f) = ~P+P- sm(9 + -9~) + ^p- cos(0~) 

From this we conclude p_ cos(#~) = p + cos(#+). 

Assume for the moment that the variables p_ cos(# - ), p + cos(# + ) / (we will rule the other 
possibilities out as minimizers below provided e is small and (3 is large). Working under this 
assumption we show that p+ must equal p_ . The sum of angles formula for sin gives 

p + sin(# + ) - p_ sin(0~) = 2e. 

as well. 

Next radial differentiation gives 

V P± = \vg - ?f cos(#+ - 9~) ± | sin(^) = 

at stationary points. Combining with the information gained from the angular differentiation, we 
find 

p T co S {8 + -9-)± € - S \n{9 ± ) = - A p ± . 
This implies that at stationary points both p + , p_ satisfy the equation 

Vg(p) ~\p = 0. 

This gives the Mean Field Equation for the standard XY model at inverse temperature (3 once all 
factors of \ have been accounted for. 

We introduce the notation ^(m) = ^g(||m||2) = <7(||m||2)— = — 3 ll 771 !!! - ^(Ilw-lbei) 
and record for reference (see [5]): 
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Proposition 4.1 Let fi c = 2. For f3 < (3 C 

d P ipfi{p) = 

has the unique solution p = where as for (3 > /3 C , there are three solutions {0, pp, — pp}. 
Further, for (3 < j3 c , p = minimizes ip while for (3 > j3 c p = is a local maximum and 
p G {p/3,— Pp} are the unique global minima. These are all the solutions to the mean field 
equation 

p = R(p) = ^ i • 

We add to this proposition the following simple observations: There exists /3q > and 5 > 
so that so that if (3 > (3q > (3 C an d Pi G {p : \p ± p^| < 5} then 

|i2(pi)-/2( / 92)|</9e-^|^i-p 2 | 
since <9 p i? = /3Varg p (f7). Further, for p ^ {p : |p ± p^l < 5}, 

Y>/3(p) - ^(p) > 4> 

and finally |p^| > ^. Assume from now on that (3 > /Jo- 
Summarizing, for all e sufficiently small (e < ^ is sufficient), f3 > (3q and if p_ cos(#~) = 

p + cos(# + ) / 0, possible stationary points satisfy p^ 1 G {Ojp^}. By assumption on (/3,e), 

local minimizers must satisfy p± = pp. We then conclude from p_ cos(#~) = p + cos(# + ) that 

9 + = ±#_ / 0. From p+ sin(6>+) - p_ sin(0 _ ) = 2e it must be that 9~ = -9+ = -9 and 

sin(#) = - p . Let us call this solution (m + , m ). 

We must still rule out the cases p_ cos(6>~) = p+ cos(# + ) = 0. It is easy to reduce to the 

two scenarios p_ = p + = or 9~ = 9 + G {§ , ^} as the other possibilities lead to larger free 

energies than these two. Both cases fall within the optimization problem 

min tpp(p). 
p 

Since (3 > /3q, we see that of the two, (m^, m^) = pp(e2,&2) has the lower free energy. Then 
the free energy difference of the latter from the absolute minimum is difference 

e 2 

0(m+,m ) - (f>(m£,m ) = - — 

by direct computation. 

Finally we consider stability of the optimizers. Again, we assume f3 > (3q. Suppose that 
£ < 5 A | and let 

= {(m + , m_) : min ||(m + , m_) ± (m + , m _ )|| < £}. 

We are interested in 4>(m,Q , ) — 4>(m + , m ) when (m^, ) ^ A^. Because we have iden- 
tified all stationary points of <p, we have that 



mm 0(niQ , m ) — 0(m + , m )= min ^H? 7 ^, m o ) ~~ <A( m+ > m )A 



r v ..« u , ...y y ^ v , , „„„ rVU' U/ T \ ' / ' 1 ey 



(m£,m ) - - {^(mJ,J?xffl + ) + <j)(R x m ,m )} = -(m^ - iJxm„ • i?x^o ~ m o )■ 
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Thus we only need to compute a lower bound on the free energy difference on dA^. 
Let Rx denote the reflection across the Si axis in M 2 . Then 

1 , ^ _ ^ .._ _ 1 
2 

This last term is positive if (mj , 771q ) G ^ and e is sufficiently small so that the Taylor expansion 
of cosine about is accurate (e < ^ will do fine). Thus to obtain a lower bound over dAg, we 
may assume m,Q = RxrriQ. Under this condition, the functional simplifies to 

<P(m+,R x m+) = X (p,0) := - -S(p) + ^ epsin0 

where p = ||7TIq" ||2 and ■ RxrriQ = cos (20). We may compute 

V 2 vf nf))= ( -l-^ p 2 S(p) + 2sin 2 psin(20)-ecos(0) \ 
XKP ' ' V psin(20) - ecos(0) p 2 cos(20) + epsin0 J ' 

The on-diagonal terms are both 0(1) while the off diagonal terms are 0(e) for f3 > /3q for 
(toq, RxtUq) € ^4g. This implies 

i?xm+) - 0(m+ m") > C{(\\m+\\ 2 - Pfi ) 2 + (9 - /3 ) 2 }. 

Since Euclidean length in polar coordinates and Cartesian coordinates are equivalent in A^ for e 
sufficiently small and (} > (3o, the proposition is proved. □ 

For m in {x : \\x\\2 < 1} the mapping m >->■ — V5(m) has an inverse, defined on R 2 . This 
inverse is given by 

v ; / e a - h du(a) 

Let 

M*(ft) = M£ e (/0 = ^ (M(f3(h + ee 2 )) + M(/3(/i - ee 2 ))) 
Note that stationary points of 4> satisfy 

m = M*(m). 

For later reference, we will need estimates on the difference 

M*(h) - M*(m) 

subject the condition \\h — m|| 2 < S. 
We have 

Proposition 4.2 77iere exists <5, eo, /?o > so that ife<eo and (3 > (3q and if 

h£{\\h' -ui\\ 2 <5}, 

then we have 

\\M*(h!) - M*(m)|| 2 < (4/3e" Cl/3 + 1 - c 2 e 2 ) \\h - m|| 2 + 0(||/ii - m|||) 
where c\, c 2 depend only on €q, (3q, 8. 



RANDOM FIELD INDUCED ORDERING 15 



Proof. Provisionally, let 5, (3q, eo be as described after Proposition 4.1. To see the stated contrac- 
tion occurs, we consider two possibilities. 
Case 1: ||m|| 2 > \\h\\ 2 . Note that 



M((3h) = R(\\h\\ 2 )h. (4.2) 

Then 

M*(h) = l - ^R(\\h + ee 2 \\)hX^2 + R{\\h-ee 2 \\ 2 )h^2) ■ 
Now by hypothesis on S 

\R(\\h ± ee 2 || 2 ) - R(\\rh ± ee 2 || 2 )| < f3e~^\\h - m|| 2 . 

Since 

R{\\m + ee 2 || 2 ) = R(\\m - ee 2 || 2 ) 
the Proposition is proved by estimating 

^\\h + ee 2 + h - ee 2 - m + ee 2 - m - ee 2 |||. (4.3) 
Letting h denote the unit vector in the direction of h and write h = ae\ + 6e 2 , we may rewrite 

\\h + ee 2 + h — ee 2 — m + ee 2 — m — ee2 1|| 

= 11^ + nT^T^ + h- Tu^ 2 -2— ei |||. (4.4) 

Let p = e 2 (||m|| 2 — ||/i|| 2 ). After a tedious perturbation theory calculation we have (4.3) 
bounded by 

||m||| [(1 - a - ap) 2 + b 2 (l- c 2 e 2 ) 2 ] + 0(||m - hg) (4.5) 

for two universal constants ci, c 2 > 0. Now /x = 0(e 2 ||rh — h\\ 2 ). Up to terms of order 0(||rh — 
hWl), ||m|| 2 ((l - a) 2 + b 2 ) = ||m - /i|| 2 so multiplying and dividing by (1 - a) 2 + b 2 we have 
(4.3) bounded by 

(1 - a - ci/i) 2 + 6 2 (1 - c 2 e 2 ) 2 _ 2 

" (i- )2 + b 2 Hl m " ^ + °(H m " ^lli)- 

By assumption on the Case, \i > 0. The required factor 1 — c^e 2 now follows whenever ||m — 
h\\ 2 < jq since a < 1 and ||rh|| 2 > \ for /3 > /3q and e < eo- 

Case 2: ||m|| 2 < ||/i|| 2 . Let /ig = pph. For /? > /?o and ||/i — m|| 2 < (5, it is always true that 

\R(\\h ± ee 2 || 2 ) - ± ee 2 || 2 )| < /3 e - c/3 ||/i - m|| 2 . 

Since \\hp — m|| 2 < — m|| 2 for ||x|| 2 = ||m|| 2 } this case reduces to Case 1. □ 

Entropy Estimates: We will need quantitative estimates on the entropy function S(m) as well 
as its finite volume approximations. Let 



( i N 

As Am) - CTGS ^£ 
I i=i 



o; - m\\ 2 < 5 
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which is a subset of Sn '■= (§1)^ and let 

Ss,N(rn) := ^logu N (As,N(m)) 

where vjj is the independent product measure on Sn spins with each one coordinate marginal 
given by Haar measure on S 1 . Both S(m) and Sg,N(m) are rotationally invariant so we may 
assume m = ph\ and consider them as functions of p, call them S(p) and Ss,n(p)- 

For each p € (—1, 1) let h(p) = argmin ft (G(/iei) — ph). Convexity implies that h(p) exists 
and is unique: 

Proposition 4.3 For each p £ (—1,1), the equation h(p) = argmin h {G(he\) — ph) has a 
unique solution which satisfies the equation 

f\ u(l - u 2 )~^e uh du 

P 



In particular, h{p) diverges as 



as \p\ t 1. 



f\{l - u 2 )-h uh du 



c - -iog(i-H) -° 



Below we will need quantitative bounds on the deviation of Ss,N(rn) from the entropy function 
S(m). The following will be sufficient. 

Lemma 4.4 Suppose that 5 > 0. Then we have 

\S(p) - Ss, N (p)\ < C5\h{p)\v 5- 2 N' 2 . 
Proof. By definition of h(p) 

e -s(,)N + s sM P)N = e °^ N ^ NMp) (AsM^ 

with Pn,H(p) the N spin probability measure tilted by h(p)e\. Because PnMp) ^ s a probability 
measure, we only need to provide a lower bound on 

PN,p(h)i A 8,N{pei))- 

The claimed bound follows from a simple application of Chebyshev's Inequality applied to A c 6 N . 

'□ 



5. Free Energy Functional Estimates for Clean Contours 

By definition, if a G X(T), the restriction of Q z (a) to each of the components 5 ex t(T), S\ n (T) 
is constant. We shall say that T is a ± contour if z (l) = ±1 on 8 ex t(T) (this notion makes 
sense by specification that T is a contour). Denoting S^ xt (T) = 5 ex t(T) in case of a ± contour 
and the empty set otherwise we let Rf (T) = J R ± n(<5* n (r)) and <5^ t (r) U (UiRf (T)) = R ± {T). 
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Given that T is a clean + contour, we say that a boundary spin configuration as(r) c i s compat- 
ible with T if 

/^rf ( e M(r)) = e r(z) for z G sp(r) and Q z (a' s{r) ) = ±1 for z G i?±(r)) / 0. 

Note here that 2 (cr^ r ^ ) implicitly takes as an argument the extended configuration (o~' s ^ , o"<5(r) c ) 
although we will continue to suppress this detail below. 
For any such cr^p^, let 

/^T) )C (°M(F)) = d r(z) for z G sp(r) and & z (a' s{T) ) = ±1 for z G i?±(r)) 

W/(r;CTAc) = ex°s(r)c 7Z , , \ ~ 

»5(r) {®z(v' s(r) ) = 1 for ze5(T)) 

where 6 x cr^p-p denotes the boundary condition with (6 x crs(r) c )z = ®z(T)o~s(r) c ,z- Let 

\\W(T;-)\\= sup W(T;a A c) 

{cr A c compatible} 

where {cr$(r) c compatible} indicates that we only consider boundary conditions for which the 
numerator does not vanish. Notions for — contours are defined similarly with + and — reversed 
and, in particular, 6 x a^c is replaced by — 6 x a\c. 
The main result which allows us to proceed is: 

Lemma 5.1 There exist 5, eo,/3o > so that if e < eo, /3 > f3 e > Po, < £ < 5 and 
p, X,k G (0, |) then the following holds: Suppose that T <Z Jj d is a (n,p)-clean contour with 
respect to to G £1 Then 

\\W{T;-)\\ < e- qiN r 

where 

qi = d {-/3(e 2 A A L- Xd )L^ d + C(/3(L A ~ 1 + ep + eL (-d/2+*0(i-A) + e -« 2 ^ A ) + 1)L( 1+A ) d } 

Let us attend to the proof of Lemma 3.3 before exposing the proof of Lemma 5. 1. 

Proof of Lemma 3.3. Let TV be fixed and consider the event X(r*, . . . , T* m , r m+ i, . . . , T m+n ) 
with sp(rj) C Aat. Then we claim 

m 

^ N (X(TI, ...,r* m , r m+1) . . . , r m+n )) < \\w(r i; .)|| 

i=i 

Once this is justified, the Lemma is proved by application of Lemma 5.1. 

The proof of this claim proceeds by induction on m. Interpreting an empty product as 1, the 
case m = there is nothing to prove, so we proceed to the induction step. Suppose the claim is 
true for any * compatible system with m = k clean contours and n dirty contours. Given a set 
{T\, . . . , r£ +1 , Tjt + 2, . . . , r^+fi+i)} of *-compatible contours and reordering as necessary, we 
may assume 

c(r fc+ i) n u*Lisp{rj = 0. 
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Assume for concreteness that is a + contour. The argument in the case of a — contour 
proceeds in a similar manner. Let {r\ , . . . , F r } denote the set of contours among {T m+ i , . . . , T m+n } 
with 5(Ti) C Int(rjfc+i) and {T[, . . . , T' n _ r } denote the rest. Then we have 

x(rj\ . . . ,r^ +1 ,r fe+2 , . . . ,r fe+n+1 ) = x(r*, . . .,r* k , r' l5 . . . ,r' n _ r ) n x(r£ +1 ,Ti, . . . ,r r ) 

Using the DLR equations we have 

A t A JV (X(ri, . . . ,T* k+l ,T k+2 , . . . ,T k+n+1 )) 

= (^(^....r^r;,...,^.^^^ (5-1) 

where {6\$= (r k+1 ){ a ) = 1} indicates that that the phase function 8 is one on 5 ex t (Xk+i) and the 
exterior ^-measurable blocks neighboring 5 e xt(Tk+i)- 
If T = (sp(r), 6 T ), let -T = (sp(r), -0 r ) and let 



r e if6 ext (Tt)c6f(r k+1 )) 

—Tt otherwise. 



Note that this transformation preserves the notion of cleanliness. It is straightforward, using the 
invariance of internal energy under spin reflections about the Y-axis, to see that 

1 i fl l^ f FHi)< ff ' El l' ll ( r Hi- f i-- f ^p(C) f irH.) 

< II W(r fc+1 ; ■ ) II 1 {<9 |,= t(rfc+i) (0 r )s i} < l X (r rfc+1 (ro,-,r rfc+1 (f .))>SK+ 1 *)t^«t(r fc+1 ) ■ t 5 " 2 ) 

Since 

1 x(r^.. ) r2JV..j"„_ji{9| i = ( ^^ 

< ^(X(H, . . • , T* k , T[, . . . , T' n _ r , T Fk+1 (fx), . . . , T Ffe+1 (f r ))) (5.3) 
the induction step is proved. □ 
Proof of Lemma 5.1: Reduction to Deterministic Weights 

For approximation purposes we introduce a block mean field Hamiltonian. In the formula B r 
denote the l < -measurable blocks in A, Bf denotes any fixed partition of B r into sets with cardi- 
nalities differing by at most one (we will call this an equal splitting below), and we let 

J e L <(z,z')^J L (r,r') 

if z G B r , z' £ B' r and B r , B' r are both i < -measurable, so that J e L K (z, z') is constant over pairs 
of l< -measurable blocks. Let us define 

mpK\<fo = -\ E jfrMtfr-SS-l E Jt^z')ai<.a e < 
z,z'eA zeA,z'eA c 



+ £ 

reAn{(2£<+l)Z} d 



zeB r 



2 ea e <> -e 2 



(5.4) 
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Let [f k denote the corresponding Gibbs measures. We introduce contour weights by 



Ai'rf (**fo(r)) = 0r{z) for z G sp(r) and @ z (a' s{r) ) = ±1 for z G fl±(r)) 



/Vr r)C (e z K (r) ) = l for ze5(T)) 



(5.5) 



defined whenever <ts(t) c is compatible with T. 



Lemma 5.2 Suppose that V is a fixed contour. For any equal splitting of Immeasurable boxes 
B r into Bp we have 



log 



W(T;a s{r) c) 



< q 2 (p,X,n)N r 



W(T; a s( r)c 

uniformly over the set of boundary conditions as(r) c compatible with T, where 

q 2 (p, A, k) = C${L X - 1 +ep + eL (-d/2+d K )(i-X)y^ 
Proof of Lemma 5.2. The key point here is that up to errors of order 

C/3 ( L A-1 +ep + eL {-d/2+dK)(l-\) )£d> 

the precise location of the local fields may be forgotten up to the information that most mea- 
surable boxes have well balanced local field statistics. 

It is easy to see that W depends on the precise choices of Bf only up to 0((3L X ~ 1 \A\). Thus 
we may work with a suitable choice of equal splitting. The choice we make is as follows: for 
each B r i < -measurable, we choose Bf- so that whenever Nf- has smaller cardinality than iV^f, 
then C Bf and then the rest of Bf is filled out by an arbitrary subset of N^. 

With this choice, we claim 

^{(ta\o- A c) - Jf(a A \o- A c) < C(X A_1 + e L( 1 - A )(- d / 2 +« d ) + e p)|A|. 

This leads immediately to the result. For future reference, we record this statement outside the 
proof along with some other internal energy approximations. □ 



For the particular equal splitting of l < -measurable boxes B r introduced in the previous proof 



let 



denote the block averages: 



*,e<,± 

T A,z 



2 

B r 



Below, it is implicitly understood that these block averages may be consider as functions on 
L°° x L°° by extending them to be t <-piecewise constant over S(T). Let be a 

two component function with each component in L°° (A) and let the boundary condition fh\c e 
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L°°(A C ). Introduce the continuum internal energy 

U^ L \rn A ,m A \fh\c) = — - j dzdz' Jl(z, z')m z ■ fh z > 

-\e j dzm+ ■ e 2 + \e [ mT z ■ e 2 - \ j dzdz',J L (z, z')m z ■ m^^ (5.6) 
1 JA 1 JA 1 JA,A C 

Proposition 5.3 Let us suppose that A is a (re, p)- clean ^-measurable subset for oj € 0. Then 

^(o-aM - ^VaKc) < c(L A - 1 + eL (i-V(-d/2+*d) + ep )| A | 
with respect to the splitting Bf introduced above. For any equal splitting, 

^Va|<tac) - U^ L \a*/<' + ,a*/<'-\a A c) < c(L x ~ l + eL (^)(-d/2+Kd) + £p) | A |_ 

Further, if (m + , m~) is an arbitrary magnetization profile in L 00 (M d ) x L°°(K d ) then 
U^{m+,ml\m A c) - U^ L \m A <' + , m A <'-\m A <)\ < cL^A]. 

Proof. This is a relatively straightforward application of the various assumptions. Let us begin 
with the claim for spin configurations: Since Jl is slowly varying, if r, r' are the centers of two 
l < -measurable boxes B r ,B' r , then 

\J L (z,z') - J L (r,r')\ < C\\ VJW^L^ 1 . 

Thus we may replace the term quadratic in the spins by 

~\ E •>t < M<fr-4 < -l E Jt(z,z')ai<.a e < + 0(L^) (5.7) 
z,z'eA zeA,z'eA c 

For any £<-piecewise constant profile, the first term is also the integral 

-\ [ dzdz'f<{z,z')ai<.ai<--\ [ dzdz* f< (z, z'V< • a[< 
1 JAxA 1 Ja,a c 

and using the approximation \ Jl(z, z') — Jl^t, r')\ < C\\V J\\l°°L x ~ 1 once again, we have 

- \ E Jl(z, z') ai< -a[<-\ E Mz, z')ai< ■ „<< = 



z,z'eA 
1 
2 



z6A,z'eA c 



\ f dzdz'J L (z,z')a e z <-a e <--l [ dzdz' J L (z, z')a e < ■ a £ < + 0(^-^1). 
z JAxA z Ja,a c 



(5.8) 



The middle claim of the Proposition follows due to the fact that the block Hamiltonian is defined 
relative to an equal splitting which is coupled appropriately to the randomness. 

For the first claim of the proposition, we still need to estimate the difference in energy con- 
tributed by the local field terms of ^{a^a^c) — Jif (ca|ca c ), whenever A is (re,p) clean, then 
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\N+\ 
N 



— 1/2 1 < i < ^/ 2+Kd except over a bad set of boxes with total number at most . Therefore 



E 

reA 

E 

reA 



E 

zeN+ 



&2 



+ Ej 



zeN r 
e 2 + 



€</<> 

2 2 



e 2 



z£B r 



■ e 2 



0(el 



-d/2+Kd 



|A|) + 0(ep|A|). (5.9) 



Combining this with (5.7) proves the first claim. Finally, the approximation for magnetization 
profiles follows a similar argument and will not be given. □ 

Proof of Lemma 5.1: Deterministic Free Energy Estimates. According to Lemma 5.2, to 
prove Lemma 5.1, it is enough for us to work with W for some arbitrary but fixed equal splitting 
of -measurable boxes B^r. From now on the notation (ai <,+ , o"f <,_ ) refers to spatial averages 
taken with respect to this equal splitting. 

Let A be a bounded ^> -measurable region in Z d . Let B K A denote the sub-sigma algebra of 

£>o,A with events determined by the block average profiles erf < ^ . An event A in Bq K a may be 
identified with a subset A* of deterministic magnetization profiles on A in an obvious way: for 



{m^,m z ) zgA 



with Hrrr^lLoo < 1, we say that (m+,m z ) G A* if (m+,m 2 



[m. 



,m. 



and this profile is taken on by (the £<-piecewise constant extension to A of) some (a z < ' , a z < ' ) 
such that a € A. Conversely, any open set A C L°° x L°° gives rise to an event A G £> \, by 
taking block averages of elements of A. Given a\c, let 



Z"(A\a A c 



/ di/(cr A )e 
J A 



-pjr"(a A \*Ac) 



For a fixed ^<-piecewise constant magnetization profile (mf,m z ), z G A and £i > 0, of partic- 
ular interest are the events 



0{m+,m z ;6) := {<r A : ||(^<' + , a%<' ) - (mj,m z ) 1 1 x,°° (A) < 

Let Dcl d bea finite union of ^> -measurable blocks, (m+, m 2 ) be a pair of vector functions 
in L°°(D) and let fhD c , z be a boundary condition on D c . We assume all functions are bounded 
by 1 in norm. Let 



Fj L ,D,eirn r ,m \m D c 



Jl(z,z') 



L 



DxD 2 

Jl(z,z') 



-m z • m 7 i 



-m r ■ rri£)c z 



/ e 2 -(m+ 
Jd 

dzS(m^) + S(m' z 



e 
2 
1 



(5.10) 



>DxD c 2 

Note that since Jl(z, z') = if \\z — z'\\2 > L, there is no loss in replacing by its restriction 
to 

3D< L := {z£D c : dist(z, D) < L}. 
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Given a spin configuration a G D n 7L d , the functional 

^J i ,D, e (^<' + ,a £ <'-|4<) 

is defined by extending the block averages from D nZ d to D in the obvious way. 

Theorem 5.4 Let A be a bounded, ^-measurable subset ofL d . Then there exists a universal 
constant c > so that 

log Z (A\a Ac ) <- (3 inf Fj ,(m + , m~ \a A i) + C{f3L x ~ l + L~3 ( - 1 ~ A - ) logL)|A| 

Also, if (rriQ z ,rriQ z ) z( za is a i^piecewise constant two component magnetization profile, and 
Ci > L~* (1 ~ A) logL 

lo g Z(0(m+ z ,m^;£i)|a A c)>^^ 

(5.11) 

Proo/ Let us prove the upper bound. To begin, we decompose {x G M 2 : ||x||2 < 1} via a 
finite collection Cl of open balls Bs{m) which intersect at most C number of times where C is a 
universal constant independent of S. Recall that h(p) = argmin fe (G(/iei) — ph) for p G (—1, 1). 
We require that this collection Cl be fine enough so that 

Sh(\\m\\ 2 ) Vt 2 ff < cL-K 1 -^ logL 

for each balls' center m. As = L 1 ~ x , 5 = L~s( 1_A ) satisfies this condition by Lemma 4.4. 
The number of balls required for is adequately bounded by cl?. 

Given A G we can find a cover C4 of A* consisting of balls Cj in L°° x L°° each with 
radius S. The centers of these balls are given by £<-piecewise constant profiles (m+ (j) , m~ (j)) 
taking values among the centers of the balls in Cl- The number of balls needed in this cover is 

2|A| 

given by, at most, \Cl\ e < ■ Since \Cl\ < cl?, this gives a total covering number of (at most) 



5|A| 

(cL) 4 



By the above and using Proposition 5.3 



\ogZ{A\a K o) < max \-/3U^(m + (j), m-{j)\a A <) + log v A (Cj) 



+ c(/3L A - 1 +Q d log L)|A|. (5.12) 

In the continuum energy, a A < is the natural £<-piecewise constant extension to A c and 

C,- = {a:(ai<'+af<'-)€C7,-}. 
Lemma 4.4 and our choice of 5 implies that 

logMCj) < cL-fM log L|A| + ^ S(m+(j)) + E S ( m * 0')) < 5 - 13 ) 
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where we used the notation A ± = U r6 A-B^- Combining (5.12) and (5.13), we have 

logZ(A\a A c) < -13 min Fj >L>A (m + (j) , m~ {j)\a%) 

+ c((3L x ~ 1 +L-t( 1 - A )logL)|A|. (5.14) 

Finally, the energy functional U is uniformly Lipschitz while the entropy function S(m) has 
the uniform modulus of continuity c||x||(l + | log ||x|||) on the domain {x : \\x\\2 < 1}, so we 
have 

logZ(A\a A c) < -p inf F JiLA {m + ,m-\a%) 

(m+,m )&A 

+ c((3L x ~ 1 +L~^ 1 - x hogL)\A\. (5.15) 
The lower bound is similar enough that its derivation omitted. □ 

Proof of Lemma 5.1: Free Energy Functional Analysis To reduce the notation, in this subsec- 
tion we do not distinguish between sets A C Z d and their extensions A C M. d . 
The following condition on ((3, e) will be important from now on: 

4/3e" Cl/3 < ye 2 . (5.16) 

where c±, C2 are the constants from Proposition 4.2. 

Lemma 5.5 Let Abe a bounded, £> measurable subset ofL d . There exist 5, eo, (3o > so that 
if {3 > (3q and < e < eo satisfy (5.16), £ < 8, p, A, k € (0, |), and L~ x < |. 77jera 



lo, 



giy(r ; (i 5( r)c) < d [-(3c(e 2 A £ 2 A L- Afl! )4 + C(/?£ A_1 + /3e~ ce2iA + L~i (1 ~ A) logL)^} 

Proo/ We proceed in two steps. Using Theorem 5.4, it is enough to obtain an appropriate upper 
bound on 

(m+,m-)eE * 

where 

= {^ (r) : ^K(r)) = M*) for z G s P( r ) and @>W = ± x for « G R ± (^)} 
in terms of 

"/ + inf , cR * jP ^' A ( m+ ' m "l c7 5( : rc)) 

(m + ,m )e_B* v ' 

for some 

B c {Q z (a' s(r) ) = l for ze 5(F)}. 
It is convenient to make a reduction. Let 



r strip = \ze 5{T) : min(dist(z,sp(r)),dist(z,,5(r) c )) > 
fstrip = | 2 € *(r) : min (dist(z, sp(r)),dist(z, ,5(r) c )) > ^| 
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All events of interest have the property that O is constant and nonzero over the connected com- 
ponents of f strip- Let 

El =E *n {(m+.mj) = @ z x (m+.m") Vz G r 5trip } . 

Under these conditions we have the following: 

Proposition 5.6 Let V be a contour. Let (f^y be any boundary condition compatible with V. 
Then there exists (vp + ,ip~) 6 E* so that 

Fj L>r (^ + ,ili-\o- e <) < inf Fj L ^{m + ,m~\a< c ) + Ce- ce2LX \F strip \. 

This is a version of a corresponding result in Chapter 4 of Presutti' s book [7]. As such, we 
postpone the proof until the end of the paper. However, from a technical perspective, this bound 
provides the main extra cost to the range of interaction (besides the intrinsic limitations discussed 
in the introduction). 

By the previous proposition, we may now work within E\ . We can bring the mean field free 
energy functional f(m + ,m~) into the picture as follows: Let 



<5(r)° = <5(r) n <^ dist(z,<5(r) c ) > 



3£> 

IT 

Because we have assumed L~ x < |, on E\ the functional F Jl ^(r) m ~ l cr f(r) c ) decouples 
into two terms: 



17 / + -I (■< \ f A A I Jl(z, z') _ 2 

Fj L ,6(r)A m > m \ a s(r)c)= / dzdz o {m z -m z >) 

1 ; J8(T)0x8(r)0 2 

+ / dzf(mt,m-)+ [ &z&z' JL ^ Z '\ m z -® z ,(T)ui) 2 
Js(r)° J5(r)Ox<5(r)o c 2 

+ C s(:r y(m + ,m~\a^ r y) (5.17) 

where 

> (m+,m-) J5(r)\5(r)° 

+ I dzdz 'Mh^l {rhz _e AT)lh f 

Js(t)\s(t)°x8(t) l 

+ / dzdz>^^(m z -fh z ,f 
Js(T)\s(r)°xS(r)\s(r) Q 1 

+ I dzdz ' J J^l {fhz - a e < T)c z , f + E(a e < r)c z ,) + 0(L^\5(T)\). (5.18) 

and the term E(a s f T y z ,) is quadratic in <7gf r y z „ independent of (m + , m~) and invariant under 
simultaneous reflections of spins about the Y-axis in each connected component of S(T) C . 
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Consider the result of reflecting m z <,+ ,m z <, ~ on <5(r)\<5(r)° and cr^c 2 on <5(T) C about 
the y-axis as necessary so that these block averages all have positive projection along the X 
axis. This transformation of magnetization profiles leaves C^pjo invariant provided we redefine 
@ x = 1 in the second term of (5.18). 

Let us extend the transformed profile to all of 5(F) by setting it to (m + , m ) on <5(T) . The 
term 

J r j L ,8(r)°,e( m+ ' rn ~) = I dzdz' JL ^ Z \ m z -m zl ) 2 

J<5(r)(»x<50(r) 2 

+ / &zf{m+,m-)+ I dzdz' JL ^ Z,Z \ fh z - 6 2 /m) 2 (5.19) 

J 6(F) ' J<5(r)°x<5(r) 0c 2 

on the right hand side of (5.17) thus represents the free energy difference between the original 
configuration and this new configuration. 

Since the new profile is in {Q z (a' s ^) = 1 for z £ 5(r)}*, we will be finished once we show: 

Proposition 5.7 With the hypotheses of Lemma 5.5, there exists a universal constant c > so 
that we have 

inf -F jL , 5( r)o, e (m + ,m-) > C(e 2 A f A L'^Nr 

(m+,m )ei?i 

Proof. Let 

Tp = {B r C sp(T) : B r is -measurable and z (T)\ Br = 0} 

= {B r C sp(T) : B r is t> -measurable, 6 z {T)\ Br = 1 

and some £>-measurable block neighbor has 6 Z (T) = —1}. (5.20) 

For B r G Tp we have t] z (T) = for some -measurable block in B r . Since the profiles 
(m + , m~) are ^< -constant, Theorem 1.1 implies 

/ dz7(m+ m-,) > ce 2 A^. 

J B r 

Otherwise, r] z (T) doesn't vanish over B r and so there are two adjacent ^-measurable blocks in 
B r for which i] z (T) takes on different nonzero values. Call them C r , C ri . As (3 > (3q and £ < |, 
we can find a universal constant c > so that 



L 



AzAz' JL{z ' Z '\ m z - m z ,) 2 > cl™L' 



jJ L (z,z') t ^ , 2 ^„ c 2d T -d 



Similarly for B r £ T ± , 

dzd/ "" v ~'~ ' {m z - m z ,y > dfL~ 
d 2 

The estimate in the statement now follows from the fact that iVp < cd(\T°\ + 1^1). □ 
This finishes the proof of Lemma 5.2 □ 
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Proof of Lemma 5.1: A Potential Flow Argument To reduce the notation, in this subsection we 
do not distinguish between sets AcZ d and their extensions A C M. d . 

Proof of Proposition 5.6. The proposition is the consequence of the following model compu- 
tation. Assume for definiteness that F is a + contour. Fix (ra^mj ) G Eq and recall that 
S ext (T) = 5(F) n Ext(T). Let 6 (T) strip = ^Strip n S ext {T). Below we restrict attention to 
S(T)stHp, but the same sort of analysis may be carried out over the components of 5i n (T). 
We may consider 

F JLMrstn P )( m i' m i\™o) 
on L°°(<5(r) strip) x L°° (5(F) strip) subject to the boundary condition induced by restricting 
(toq, nriQ ) to 5(F) n 5(F) c Strip and subject to the constraint that the profile 



(m + ,m ) : = 



(m\,m l ) on S(T)stri P , 
(mQ,mQ ) otherwise, 



lies in E'q. By definition, this is equivalent to requiring that @(mf , m^) = 1 on 5(F) strip- 

For each m2 G {x G M 2 : ||x||2 < 1}, let us define M _1 (7712) = — VS(m2). Because S(m2) 
is strictly concave in this region, the inverse M(h) is well defined. It takes the explicit form 

*<») _ 

v ; / e (jh Av(a) 

for all liel 2 . 

We consider the differential equation: 

$m± = - [m* - M(0(J L * in ± ee 2 ))] (5.21) 

on <5(T) s tr jp subject to the boundary conditions 

(m+(t),mj(t)) = (m+m -) for (t, z) G {t = 0} x <5(r) 5trip U {t > 0} x 5(F) c strip . (5.22) 

Here, Jl * fa is the convolution of with m: 

JL*fh z := I JL(z,r)fh r dr. 

JR d 

The stationary solutions to this differential equation satisfy 

m± -M(/3(J L *m z ±ee 2 )). (5.23) 
Note that (m + , m + ) solves the vector equation 

m ± = M(/3(m±ee 2 )) 
so we may think of (5.23) as a kind of generalized mean field equation. 

Lemma 5.8 There exist 5, eo, fio > so that if (3 > (3q and e < eo satisfy (5.16), £ G (0, <5) a«<i 
(j?Iq ,mg) G Eq. Then we have existence and uniqueness of (5.21), (5.22) for all time t G [0, 00) 
in L°° (6(F) strip) x L°° (5(F) strip)- Further, for each initial condition (m^mg) G Eq, there 
is a unique stationary solution to (5.23) with boundary condition (tHq, tHq )\s(r)° st ■ I ts £<- 
measurable block average profile lies in Eq. 
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Proof. For a start let us assume that 5, eo, Po > as in Proposition 4.2. We will adjust (3q in a 
moment. The function m >->■ M(f3(m ± €62)) is C-Lipschitz within {x : ||x — m||2 < £} since 

V ft M(/?/i = pVM(ph) = ((a-a) - (a)-(a)) 

where : denotes the outer product of two vectors and the expectation is taken under v tilted by 
e /3cr /i j n f act ^ mcreasm g ^ Q as necessary 

||/3VM(/3J*m±ee 2 )|| < 1 
on E since in this case (3 J * m ± ee2 are uniformly bounded away from 0. So if 

K ± (m + ,m~) = m ± - M((3( J * m ± ee 2 )) 

then if ± are (at worst) 2-Lipschitz (in L°°). We will use this at the end of the proof as well. 

Using the Lipschitz continuity, we have global existence and uniqueness on L°° x L°° for any 
initial condition in Eq . Let 

E = {(m + ,m~) € L°° x L°° : \\fh - m|| L oo < £, ||m ± || L oo < l,m ± • e x > 0}. 

Note that Eq C E. We first claim that the dynamics preserves E: 

By continuity in t, given To sufficiently small and positive, we may assume that on the interval 
[0, To] the solution is in E. Note that the constant profile (a + , a~) = (m + , m~) is in Eq. Since 
this profile is stationary we have 

^d t \\d(t) - m 2 (0ll2 = — 1 1 «• - ^2 ||i + [« - rhi] ■ 

-(M{f3{J L * 6 + ee 2 ) + M(/3( J L ★ a - ee 2 )) - -{M((3(J L *fh 2 + ee 2 )) + M(/3( J L * m 2 - ee 2 ))) . 

(5.24) 

By assumption on /3 and Proposition 4.2 

||^(M(/3(J L *a+ee 2 )+M(/3(J L *a-ee 2 ))^ 

< (l-ce 2 )||J L *a-J L ★m 2 || 2 (5.25) 

for some c > 0. Since J* is a contraction in 

^<9t||a - m 2 \\l = -\\a- m 2 ||| + (1 - ce 2 )£||a - m 2 || 2 

for all t < T). By a continuity argument, if a priori ||a(0) — m 2 (0)|| 2 < ^ it remains so for all 
times t G [0, 00). 

Consider next the boundedness of the individual components m^. The condition mf -Si > is 
evidently preserved by considering the flow velocity along with the fact that ||m— d\\2 < > 0. 
If the initial condition satisfies ||m^||£,oo < 1, then since \\M(/3( Jl *rri2 ± ee 2 )))|| 2 is uniformly 
bounded away from 1 (on E), it is easy to see that the dynamics is a coordinate wise contraction 
when either component lies in p < \\rri^ || < 1 for some p sufficiently close to one. We conclude 
via another continuity argument that the dynamics preserves E. 
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Under the flow, the free energy functional Fj L ^^ strip (m + (t),m (t)\mo) subject to initial 
conditions lying in E, evolves as: 

d tFj L ,S(T) strip = 

dr[j-km + ee 2 + VS(m + )] ■ [m + - M(f3( J -kfh + ee 2 ))] 

+ \ j dr[J -km - ee 2 + VS(m~)] ■ [m~ - M(/3( J -km- ee 2 ))] (5.26) 

Now M and — VS are inverses of one another when the latter is restricted to ||m|| 2 < 1, so each 
integrand is of the form 

-(M(m) - M{h)) ■ (m-h). 
Since VM > as an operator, we conclude that dtF.j L ,s(Y) st ri P * s decreasing. It is strictly 
decreasing unless (m + , m~) satisfies (5.23) because 

v ■ VMv = Vw h {a -v)^0 

for any unit vector v and any li£l 2 . Here Var^ denote the variance of a random variable under 
the measure defined by tilting v by e a h . It is easy to see that the functional is bounded below 
on E, so we conclude that any limit point of the flow started from an element of E must be a 
stationary point. 

From (5.23), we have 

m ± = M{P(J L *m± ee 2 )) 

so that the two components of any stationary point MUST be determined by the average vector 
m. In particular using the representation (4.2) 

\\mf - m ± || 2 = \\M{f3(J L -km z ± ee 2 )) - M(/3(J L ★ m ± ee 2 ))|| 2 < \\m z - m|| L oo 

where we have used that M(j3v ± ee 2 ) is 1-Lipschitz in v and that J^* is a contraction in L°°. 
We have proved that any stationary solution has block averages in Eq . □ 

The next order of business is to show that if mf is a stationary point with boundary condition 
m$ then m r is very close to m if r is deep inside 5 (T) strip- 

Lemma 5.9 With the same hypotheses as in the previous proposition, let {mf , m~ ) be a station- 
ary solution to (5.23) on 5(F) strip subject to the boundary condition coinciding with (m^, m^) € 
.Eo on 5(F) strip- Then there exist universal constants C,c> so that 

\\fh r — m r || 2 < Ce~ ct2LX 

whenever d(r, S(F) c Strip ) > l -f 

Proof. If m G E is stationary then (using notation from Section 4), 

fh - m = M*p(m) - Mg(m). 

Hence by Proposition 4.2, if on some region D \\m — ra\\ L oo^ D ) < £ then on the interior D° L = 

{x G D : dist(x,D c ) > L} 

\\m - mlloo < (1 - ce 2 )£. 
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Iterating, if dist(r, D c ) > kL then 

\\m r - m|| 2 < (1 - ce 2 ) fe £. 
Hence we conclude that any stationary solution has 

\\fh r - m|| 2 < Ce" ce2LA ^ 

if dist(r, 5(r)g tri ) > -g 1 . Since the operator norm of VMg is bounded by 1 the same holds for 
rrv^ using (5.23). □ 

Finally, if (m+,mj) is a stationary solution to (5.23) on 5(T)stHp subject to the boundary 
condition coinciding with (m}, mg ) G i?o on <$(r)g tri , then its block averages (m^'" 1 ", m^' - ) 
lie in Eq and satisfy 

||mf<' ± -m ± || 2 <Ce- ce2LA 
if dist(r, $ (T) strip) > F rom tne proof of Lemma 5.8 

F ^,a(r) s ^ P ( m+ > m ~l™o)= inf F Jii5(r)5trip (m + ,m-|7fi ) 

(mj",mj )&E 

Note that 

|| J -kfh z — J *m<\[2 < cL A_1 

and that 

/ dz/(m+ mj) > e d < f(m t <' + ,m e <'-) 

J\z-r\<e< 

on E'q since / is convex on {||(m + , m~) ± (m + , m~)|| < £}. 
Thus we have 

J? Ji, ) «(r)stri J ,( m<< ' + ' m/< '"l"»o) " F ^,5(r) st , iP ( m+ ' m_ No) < cL A - 1 |<5(r) 5 trip|. 

To summarize, we have shown the existence of an l < -constant magnetization profile (m 2 < ' + , 
so that 

< inf F Jii<5(r)strip (m + ,m-|mo) + cL A - 1 logL| ( 5(r)^„ P | (5.27) 

(mi.mj )eE 



and 

||mf<' ± -m ± || 2 <Ce- ce2LA 



if dist(z, S(T)g trip ) > -g 1 . To complete the argument, we simply modify m z < ' on &\st(z,5(T) c Strip 
^ to be equal to m^. Call the new profile (mif 2 ' + , mif 2 '~). Then 

f Jl ,«(r) Strip (m£< •+ , m£< | m ) 

< inf F JlMT) (mt,mi\mo) + cL x ~ l \8{T) stnp \ + Ce^^ \5{T) stnp \ (5.28) 

(m{,m 1 )eE 

since the on-site term in Fj L ,s(r) strip can on ly De decreased by this modification. This concludes 
the proof of Proposition 5.6. □ 
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